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Abstract 

We formulate the single transverse spin asymmetry in heavy quarkonium production in lepton- 
nucleon and nucleon-nucleon collisions in the non-relativistic limit. We find that the asymmetry is 
very sensitive to the production mechanism. The final state interactions with the heavy quark and 
antiquark cancel out among themselves when the pair are produced in a color-single configuration, 
or cancel out with the initial state interaction in pp scattering when they are in color-octet. As a 
consequence, the asymmetry is nonzero in ep collisions only in the color-octet model, whereas in 
pp collisions only in the color-singlet model. 
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1. Introduction. Single-transverse spin asymmetry (SSA) is a novel phenomena in 
hadronic reactions , and has long been observed in various processes. In these pro- 
cesses, a transversely polarized nucleon scatters off a unpolarized nucleon (or virtual photon) 
target, and the final observed hadron have asymmetric distribution in the transverse plane 
perpendicular to the beam direction depending on the polarization vector of the scattering 
nucleon. Recent experimental studies have motivated much theoretical developments for 
understanding the underlying physics associated with the SSA phenomena. 

An important theoretical progress made m the last few years is the uncover of the crucial 
role of the initial/final state interactions [3], which leads to a non-vanishing SSA in the 
Bjorken limit in the semi-inclusive hadron production in deep inelastic scattering (SIDIS) 
and Drell-Yan lepton pair production processes. These initial/final state interactions were 
understood as a result of the gauge link in the gauge invariant transverse momentum de- 
pendent (TMD) quark distributions , and the SSA can be traced back to a naive 
time-reversal-odd distribution, the so-called Sivers function [3]. The difference between the 
initial and final state interaction in the above two processes results into a sign change for 
the associated quark Sivers function and the SSAs. The Sivers function has been extended 
to the gluon sector, and the associated prospects have been investigated in recent years 
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H eavy quark and quarkonium productions are natural probes for the gluon Sivers func- 



tion 
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14| . Especially at low transverse momentum, their production will be sensitive 



to the intrinsic transverse momentum [15[, and can be used to study the gluon Sivers ef- 
fect. Meanwhile, the heavy quarkonium production has attracted many experimental and 
theoretical investigations in the last decade, starting from the anomalous production dis- 
covered at the Tevatton pp expeti„.ent Q and a theoretical f_tk tj. =tudyt„g the 
heavy quarkonium system, the so-called non-relativistic QCD (NRQCD) [17|. The basic 
argument for the NRQCD factorization is that the heavy quark pair are produced at short 
distance in a color-singlet or color-octet configurations. The hadronization of the pair (in 
singlet or octet) is described by the associated matrix elements, which can be character- 
ized according to the velocity expansion [17]. This framework had success in explaining 
some experimental observations. Yet, we do not have a definitive answer for the production 



mechanism 
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In this paper, we formulate the SSA in heavy quarkonium production by carefully exam- 
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ining the initial and final state interaction effects. We follow the general arguments of the 
NRQCD factorization. The SSA depends on the color configuration of the pair produced 
at short distance, although the hadronization details are not relevant to obtaining a non- 
zero SSA. Therefore, the experimental study of the SSAs shall shed light on the production 
mechanism for the heavy quarkonium. We will focus on the gluon-gluon (photon-gluon) 
fusion contributions in two processes: one is the lepton-nucleon ep scattering and one is 
nucleon-nucleon pp scattering, schematically as 

a{Pa, s^) + b{Pb) ^ [gg](^'') + x^h + x , (i) 

where A is the polarized nucleon with momentum P4, B the un-polarized hadron or photon 
with momentum Pb- The heavy quark pair QQ produced at short distance (~ 1/Mq, Mq 
the heavy quark mass) can be in a color-singlet (labeled with (1)) or color-octet (labeled 
with (8)) configurations, and H represents the final physical quarkonium state. The SSAs 
coming from the initial/final state interactions are calculated based on the interactions 
between the quark pair or the incident gluon with the remanet of the polarized target. We 
are interested in a kinematic region of low transverse momentum for the heavy quarkonium, 
P± <^ Mq, where the TMD parton distribution is relevant. In the following, we will present 
a general analysis for the SSAs coming from these interactions in this kinematic region. A 
detailed argument for low P± factorization in terms of the TMD parton distributions is also 
important. Especially, when P± is the same order of the heavy quark kinetic energy in the 
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quarkonium's rest frame, the hadronization of the heavy quark pair might inter: 
initial gluon radiation and break the TMD as well as NRQCD factorization 
will come to this important issue in the future. 

2. Final state interactions with the quark pair in the non-relativistic limit. We 

first discuss how to formulate the final state interactions effects with the quark pair, taking 
the lepton-nucleon scattering as an example. In this process, the dominant subprocess is 
the photon-gluon fusion channel, where the photon comes from the lepton radiation. As 
we discussed in the introduction, in order to obtain a nonzero Sivers-type SSA, we have 
to generate a phase from either the initial or final state interactions. Because the lepton 
(photon) is colorless, there is no initial state interaction in this process. 

In Figs. 1 and 2 we plot the generic final state interaction diagrams, for color-singlet and 
color-octet cases, respectively. If the pair are produced in a color-single configuration, a 
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FIG. 1: Vanishing SSA in photo (lepto) -production of heavy quarkonium when the heavy quark pair 
are produced in a color-singlet configuration. The two final state interactions with the quark (a) 
and antiquark (h) cancel out each other; the final state interactions with the unobserved particles 
cancel out among different cuts. 



gluon has to be radiated. In general, we will have final state interactions with the quark 
pair (Fig. l(a,b)) and the radiated gluon (Fig. l(c,d)). The interaction with the unobserved 
particle (here is the radiated gluon) vanishes after we summing different cut diagrams [2], [2^ . 
For example, the contribution from Fig. 1(c) is proportional to pole contribution from the 
propagator labeled by a short bar with momentum k' -\- ki, 

1 



(k' + k^y + ie 



6{{k'f)\^ole = -^7r6{{k' + k^Y)6{{k 



f\2\ 



whereas the contribution from Fig. 1(d) will be 

1 



{k') 



m + A;i)')|poic = +i7i5{{yf)5{{k' + hf) 



Clearly, these two contributions cancel out each other, because the other parts of the scat- 
tering amplitudes for these two diagrams are the same except the above pole contributions 
which are opposite to each other. The above result is quite general, and shows that all the 
final state interactions with the unobserved particles do not contribute to the SSA for the 
associated process. 
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Therefore, we only need to consider the final state interactions with the quark pair. Two 
example diagrams are shown in Figs. 1(a) and (b). In the non-relativistic limit, the final 
state interaction with the quark in Fig. 1(a) can be derived as follows. 



9 ^n. P.,,^rT.b 



^k--;ij)T''r , (2) 



-kt +ie ' 2 

where k and P are momenta for the quark and the quark pair, respectively, $ represents 
the wave function for the pair, ij are color indices for the quark and antiquark, b is the 
color index for the gluon attaching to the quark line, p is the index contracted with the 
gauge potential A'^, and F represents other hard part for the scattering amplitude. The 
light-cone momentum components are defined as /c^ = {k^ ± k^) / and we assume that 
the polarized nucleon is moving along the +z direction: Pa = (P4,0~,0_l). In the non- 
relativistic limit, each of the quark and antiquark carries half of the pair's momentum, 
k ~ P/2, and they are on mass shell: k'^ ~ Mq. More over, the dominant contribution from 
the gluon interaction with the nucleon remanet comes from the gauge potential A'^ in the 
covariant gauge, and where the gluon momentum is coUinear to the polarized nucleon. In 
above derivation, we have also only taken the leading power contributions, and neglected all 
higher order correction in terms of 1/Mq. This derivation shows that we can simplify the 
final state interaction as an eikonal propagator. Similarly, when the gluon attaches to the 
antiquark. Fig. 1(b), the contribution will be proportional to 

-9 _T''^k-^;zj)T , (3) 



-k+ + ie ' 2 

where the F is the same as above. The minus sign of this result comes from the interaction 
with the antiquark. We can combine the contributions from these two diagrams together, 
and it is proportional to 

9 - ij)T' - T%{k - V . (4) 



-kt + ie\''" 2 

From this, we immediate find that the SSA contribution vanishes if the pair are produced 
in a color-singlet configuration, for which the color matrix for $ will be oc 5^, and the 
above two terms cancel out each other completely. This observation is generic, and only 
depends on the non-relativistic approximation we made and it is valid for higher orders 
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FIG. 2: Nonzero SSA in ep collisions arise when the heavy quark pair are produced in a color- octet 
configuration. The two final state interactions with the quark (a) and antiquark (h) have net effects. 

too. Therefore, we conclude that the SSA for heavy quarkonium production in ep scattering 
vanishes in the non-relativistic hmit in the color-singlet model. 

However, the above arguments do not hold when the pair are produced in a color-octet 
configuration, for which we will have net effects from the two final state interactions. We 
show these interactions in Fig. 2. Because the pair are in color-octet, the color matrix of 



the wave function can be parameterized as $ 



(8,c) 



oc Tf, where c 



, 8 represents the 



color-index of the pair. The sum of the two final state interactions will be proportional to 

—^{~^hc,)^^''''^T , (5) 
—k{ + le 

where the latter factor is the same hard part in the above scattering amplitude. This result 
can be summarized into a diagram shown in Fig. 3(a), and can be easily extended to a 
two-gluon exchange contributions (b), which is proportional to 

g 9 



lfdce){-lhef)^^'^^^T . 



(6) 



— ^2 + — k2 + ie 

When it is generalized to all orders, we will find the final state interactions can be summed 
into a gauge link associated with the gluon distribution, and the SSA depends on the gluon 
Sivers function, which is the spin-dependent part of the following distribution 25| . 



-ixP+^-+ik_L-ii 



P+(27r)3 

x{PS^\F+^{r,U)4~£ , -^0,0, F^+{0)\ PS, 



(7) 



where sum over color indices is implicit. F^'^ is the gluon field strength tensor, F^^'^ = 
d^A'^ — A'^ — gfabcA'^A'^, x the momentum fraction carried by the gluon, k± the transverse 
momentum. the process-dependent gauge link 0, For the diagrams in Fig. 3, the 
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FIG. 3: Summarize the two final state interactions in Fig. 2 into a gauge link associated with the 
gluon interaction. 

gauge link sums all the final state interactions with the quark pair (in color-octet state), for 
which we have future pointing gauge link going to +00, 



Pexp 



(8) 



in the covariant gauge, where = A^f^ is the gluon potential in the adjoint representation. 



with t 



ab 



—ifabc- A transverse gauge link at the spatial infinity is needed to retain the 



gauge invariance in a singular gauge 
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A potential contribution to the above gluon Sivers function comes from the quark Sivers 
function by splitting. Following the calculations in I24I] , we find that the large k± gluon Sivers 
function can be generated from the twist-three quark-gluon correlation function Tp{x) j^. 
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as e^^SlklN, 



IT 



x' \\ dx' ^ ' ' 



X 



2vr2 (A;2)2 2 

1-0 [1 + {1-0']+Tf{x',x')2{1-0'} 



(9) 



where ^ = x/x' and a sum over all quark flavor is implicit. The overall sign of this distribution 
depends on the relative contribution from the derivative and non-derivative terms, and the 
up and down quark contributions. 

In summary, the SSA in ep collisions vanishes in the color-single model, but survives in 
the color-octet model. 

3. SSA in pp collisions. Now, we turn to the SSA in heavy quarkonium production in 
pp collisions. In this process, we have both initial and final state interactions. However, as 
we showed above, when the heavy quark pair are produced in a color-singlet configuration, 
there is no final state interaction, and we only have initial state interaction contribution. 
We show a typical diagram in Fig. 4(a) in the gluon-gluon fusion subprocess. This diagram 






FIG. 4: Only initial state interactions contribute to the SSA in hadron-production process in the 
color-singlet model (a). On the other hand, the SSA vanishes in the color-octet model, because the 
cancelation between initial (b) and final (c) state interactions. 



is in particular the dominant channel for Xc production in the color-singlet model 17|]. The 
initial state interactions for this diagram can be analyzed, and we find the SSA will depend 
on the gluon Sivers function defined in Eq. ([7]), however with a gauge link going to — oo, 



r'—oo 



Pexp 



(10) 



and the associated SSA will be opposite. 

However, if the quark pair are produced in a color-octet configuration, we will have 
both initial and final state interactions. After taking into account both contributions, the 
associated gauge link becomes, C^-^^^^ — > which is responsible for the SSA. However, 

the contributions from the two gauge links cancel out each other completely. This is because 
the combined gauge link is invariant under time-reversal transformation: the two gauge links 
transform into each other and the combined one remains the same. We have also checked 
this cancelation by an explicit calculation up to two gluon exchange contributions. This 
result indicates that a standard TMD factorization breaks down for this case, similar to 
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recent discussions on the dijet-correlation in hadronic collisions 

4. Summary. In this paper, we have formulated the single spin asymmetry in heavy 
quarkonium production in high energy scattering in the non-relativistic limit which should 
be valid in the limit Mq —>■ oo. Very interesting observations were found. The SSA vanishes 
in ep collision when the pair are produced in a color-singlet configuration, and a nonzero 
SSA arises from the color-octet contribution from the gluon Sivers function with the gauge 
link pointing to +oo. On the other hand, the SSA in pp collisions in the color-single model 
depends on the initial state interactions leading to a gluon Sivers function with gauge hnk 



pointing to — oo, and the initial and final state interactions cancel out each other in the 
color-octet model and the SSA vanishes to all orders. 

In our discussions, we only considered the gluon-fusion contributions. The quark channel 
is also important for heavy quarkonium production, especially in pp collisions at relative low 
energies. In this case, the SSA will depend on the quark Sivers function from the initial 
state interaction in the color-single model, which however is opposite to the Drell-Yan SSA 
because of different color-factor, with a relative factor {—1/2Nc)/Cf- In the color-octet 
model, both initial and final state interactions contribute. At one-gluon exchange order, the 
final state interaction dominate over the initial one, and their total contribution carries an 
overall factor —{Nc/2 + 1/2Nc)/Cf relative to that in the Drell-Yan process for each quark 
flavor contribution. It will be interesting to compare the SSA in these two processes, and 
scan the energy dependence 
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